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The performance of magnetic nanoparticles is intimately entwined with their structure, mean size and 
magnetic anisotropy. Besides, ensembles offer a unique way of engineering the magnetic response by 
modifying the strength of the dipolar interactions between particles. Here we report on an experimental and 
theoretical analysis of magnetic hyperthermia, a rapidly developing technique in medical research and 
oncology. Experimentally, we demonstrate that single-domain cubic iron oxide particles resembling 
bacterial magnetosomes have superior magnetic heating efficiency compared to spherical particles of 
similar sizes. Monte Carlo simulations at the atomic level corroborate the larger anisotropy of the cubic 
particles in comparison with the spherical ones, thus evidencing the beneficial role of surface anisotropy in 
the improved heating power. Moreover we establish a quantitative link between the particle assembling, the 
interactions and the heating properties. This knowledge opens new perspectives for improved 
hyperthermia, an alternative to conventional cancer therapies. 



Consider an assembly of single- domain particles with uniaxial anisotropy, coupled to each other by the 
magnetic dipole-dipole interaction. Such interactions can be tuned by adjusting the size, the magnetiza- 
tion, and the volume fraction of the particles \ Even in the superparamagnetic regime, the collective 
magnetic behavior will differ from that of isolated particles^. A potentially interesting area where this effect 
may find applications is biomedicine. An example of the latter is the optimized design for multicore particles 
achieving enhanced transverse relaxivities for magnetic resonance imaging^'^. Yet magnetic nanoparticle suspen- 
sions have gained an important role in cancer treatment with AC hyperthermia. In synergy with chemotherapy or 
radiotherapy, selective targeting and localized heating of tumor cells can be tuned, leading to modalities with 
shorter time regimes even with a lower dosage^. Indeed, the efficiency of this type of radio -frequency thermo- 
therapy has been demonstrated on several types of cancers including brain tumor, prostate cancer, and invasive 
breast carcinoma^. Although encouraging results on palliative care indicate that even non-optimized particles 
with the appropriate size distribution can deliver adequate heating power if present in sufficiently high concen- 
trations^, concerns have been raised regarding the toxicity for cancer-directed therapy^. In order to minimize the 
potential side effects arising during the clinical treatments, the quantity of nanoparticles administered needs to be 
as small as possible but still retaining the desired effect (see Supplementary Information). For this purpose, to 
reach the therapeutic temperature with minimal particle concentration in tissue, the magnetic nanoparticles 
should exhibit high inductive specific absorption rate (SAR). 

This quantity depends on the nanoparticles' properties, such as mean size, saturation magnetization (Ms) and 
magnetic anisotropy (K), but also on the alternating magnetic field amplitude (Hmax) and frequency (/). In 
previous work^, heating has been predicted for superparamagnetic nanoparticles within a model in which SAR 
primarily depends on magnetic spin relaxation processes. It was shown that the crossover between Neel and 
Brown regimes of relaxation depends on the anisotropy constant and particle volume^, thus defining, for each 
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frequency, a narrow range of K and size values for optimal SAR. 
Furthermore, note that the aforementioned model is suitable to cal- 
culate SAR only for very small nanoparticles in the diluted regime at 
low magnetic fields whereas for clinical applications the relative 
influence among particles cannot be neglected (different results indi- 
cate that the nanoparticles are generally agglomerated in in vitro and 
in vivo experiments). For instance, results for dense systems have 
shown that dipolar interactions not only affect the susceptibility^ 
but also the blocking temperature transition^^, and the motion of 
particles in solution Hence, the selection of the most advantageous 
materials for clinical hyperthermia treatment is still a matter of 
debated 

Previously, the tuning of the magnetocrystalline anisotropy and 
its influence on magnetic heating efficiency has been issued by Lee 
et aV^ in exchange- coupled nanoparticles. Note that K can also be 
controlled by changing the nanoparticle shape. Yet, another inter- 
esting option to improve the heat generation from magnetic nano- 
particles could be to explore single domain particles displaying 
hysteresis losses^^, especially if one considers the SAR reduction usu- 
ally found in the small superparamagnetic particles after immobiliza- 
tion (for instance, into cells) Magnetization reversal calculations 
within the Stoner-Wohlfarth model for independent, randomly 
oriented uniaxial single- domain particles, lead to losses scaling 
approximately with MsHc, being He their coercive field. In this 
regard, metalUc iron-based nanoparticles (instead of iron-oxides) 
are the best candidates as they virtually present the highest saturation 
magnetization. For instance, iron nanocubes^^ (with effective aniso- 
tropy constant Keff = 9.1 X 10^ J/m^ and Ms = 1.7 X 10^ A/m) 
display the highest specific losses reported in the literature so far, 
provided that the field amplitude employed is sufficiently high to 
ensure the remagnetization of the particle (SAR about 3000 W/g at 
l^oHmax ^ 73 mT and / = 274 kHz). Nevertheless, soft magnetic 
nanoparticles seem to be preferable for the purpose of hyperthermia 
within the range of magnetic fields used in clinical trials^^. On the 
other hand, the low chemical stability at physiological conditions and 
adequate biocompatibility are drawbacks of metallic particles, thus 
forcing the combinatory use of complex core -shell architectures^°'^\ 
Consequently, up to date, the majority of investigations focused on 
magnetic iron oxides Fe304 (magnetite) and y-Fe203 (maghemite), 
which have been proved to be well tolerated by the human body. 
Therefore, it is highly desirable to determine hyperthermia efficiency 
of iron oxide particles with mean size ranging above the superpar- 
amagnetic limit (> 15 nm) and below the optimal size (<50 nm) for 
internalization into mammalian cells^^, especially if one has in mind 
the case of intracellular hyperthermia^^. 

Inspired by the excellent heating properties of bacterial magneto- 
somes having a mean diameter of the magnetite crystals of about 
30 nm^^'^^, here we demonstrate that cubic nanoparticles possess 
superior heat induction power compared to spherical iron oxide 
particles of similar size. This higher hyperthermia performance is 
mainly related to the higher surface magnetic anisotropy and the 
tendency towards aggregation into chains facilitated by the cubic 
shape. Although previous experimental studies of the magnetic pro- 
perties of cubic particles have been done^^'^^, showing that cuboids 
could be good materials for the biomedical imaging applications^^'^^, 
little has been done to compare the effects of spherical versus cubic 
particles in hyperthermia^^. Some of us have previously reported 
dependence of the SAR on the mean particle size^°'^\ Based in our 
know-how in the size and shape control of iron oxides, here we 
provide a fundamental aspect of tuning the heat dissipation effi- 
ciency by adjusting the morphological profile of magnetic particles. 
For our proof- of- concept, we first synthesized, structurally and mag- 
netically examined 20 nm cubic particles as prototype hyperthermia 
carriers in contrast to standard spherical ones. We have also com- 
pared the 20 nm cubic particles to cuboids with an increased size of 
40 nm in order to account for thermal contributions to the energy 



barrier for magnetization reversal. Theses sizes were deliberately 
selected so that the particles safely reside within single magnetic 
domain range^^'^^. On the one hand, by performing systematic AC 
magnetic measurements in aqueous solution, we demonstrate that, 
the SAR is strongly related not only to the shape, the volume and the 
concentration of particles, but it might as well be influenced by their 
aggregation into chains. On the other hand, the analysis of the 
numerical calculated hysteresis loops allows for the correlation of 
their magneto -structural properties, the influence of anisotropies, 
and dipolar interactions, with hyperthermia measurements. 

Results 

Figure 1 shows the transmission electron microscopy (TEM) view of 
representative aggregates of cube-shaped iron oxide particles (see 
Supplementary Fig. SI for a TEM image of the spherical particles). 
The difference in contrast within the image is because the different 
crystallographic orientations of individual particles with respect to 
the electron beam. In this regard, there is no apparent contrast vari- 
ation within each nanoparticle, thus pointing to the fact that particles 
are completely oxidized during synthesis A high- resolution 
(HRTEM) image clearly attests to the monocrystalline structure with 
a lattice fringe of 0.24 nm characteristic of (222) planes of the spinel. 
As can be seen in Fig. lb (inset), the cubes organize themselves into 
chains with sharing flat faces of the {100} type. 

We used electron tomography to generate a 3D illustration of the 
self-assembly of nanoparticles (see Fig. 2), similar to that described in 
biogenic magnetosomes. The faceted cubic particles are recon- 
structed separately in Fig. 2c (see Supplementary Movie SI and 
Fig. S2 onUne for details). It should be noted here, that van der 
Waals attractions are expected to be weak in our nanoparticle solu- 
tions due to the thick surfactant layers (the decanoic acid 
CH3(CH2)8COOH has a carbon chain length of 1.4 nm). The ther- 
mogravimetric analysis (TGA) was employed to verify the coating 
structure of the nanoparticles. The experimental values of weight 
losses (see Supplementay Fig. S4) are in reasonable agreement with 
the theoretical estimation assuming that the surfactant forms a close- 
packed monolayer on the nanoparticles^^. Thus, the formation of 
such chain-like aggregates is favored because the energy of the mag- 
netic dipole-dipole interaction is presumably larger than the thermal 
energy, and the van de Waals or electrostatic interactions within the 
solution^^. 

A simple way of visualizing the consequences of these dipole inter- 
actions is to look at the magnetic response. Field and temperature- 
dependent magnetic measurements were performed on the dried 
crystalline powder using a superconducting quantum interference 
device (SQUID) magnetometer. Noteworthy, the remanent magnet- 
ization (Mr) values (see Fig. 3, and Supplementary Fig. S5 for hys- 
teresis loops of the spherical particles) are well below the 0.5 Ms 
expected for magnetically independent uniaxial Stoner-Wolhfarth 
nanoparticles, thus signifying non-negligible dipolar interactions 
between the particles. The 40 nm cubic particles are ferrimagnetic 
at room temperature with a saturation magnetization value up to 
89 Am^/Kg, and coercive field of 5.5 mT. The smaller the particles, 
the smaller the magnetization of the sample is, probably because of 
the appearance of cation vacancies and surface spin canting on 
decreasing the particle size^^"^°. Hysteresis loops also indicate a 
reduction in coercivity as the mean particle size is decreased, and 
an increase of both Ms and He as the temperature is reduced from 
300 K (Fig. 3b). Importantly, data suggest that the 20 nm particles at 
room temperature are in the transition regime between superpara- 
magnetism and a blocked state. Both observations are consistent with 
a lowering of the energy barrier for magnetization reversal that leads 
to faster relaxation by thermal fluctuations^^ The effect of temper- 
ature upon the magnetic anisotropy is a topic for over 50 years^^. 
We have evaluated the effective anisotropy constant from the 
law of approach to magnetic saturation (see the Supplementary 
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Figure 1 | TEM images. Iron-oxide nanocubes (a) with average edge 
length of 20 ± 4 nm; inset reveals 2D self-assembly arrangements, 
(b) Corresponding TEM micrograph of 40 nm nanocubes. As can be seen 
in the larger area view, the particles organize themselves in different 
chain-like configurations, (c) High-resolution observation of crystal 
structure revealing (222) fringes of the inverse spinel iron oxide. Inset 
exhibits the EFT spectrum. 

Information for details)^^ The parameters of these particles are sum- 
marized in Table SI. Remarkably, the Kgff can be modulated by 
varying the size and shape of the nanocrystals, with values above 
those for bulk Fe304 (about 1 1 kj/m^) or y-Fe203 (13 kj/m^) at room 
temperature. We deem cubic particles exhibit higher anisotropy 
energy values than that for spheres due to shape contribution. 




Figure 2 | Tomography. 3D reconstruction of cuboctahedral shape 
particles, from images obtained at different tilt angles relative to the 
electron beam, after 40 iterations (see Methods), (a) Nanocube cluster. 
Neighboring nanocubes have {100} surfaces face to face separated by a 
distance above 2 nm due to hydrocarbon ligands. (b) Single nanocube in 
its original context, (c) Illustration of small deviations from perfect cubic 
symmetry. 

Note that the sphere has the smallest surface area among all surfaces 
enclosing a given volume. It is therefore not surprising that we find an 
increased anisotropy in the case of cubic particles, compared to the 
spherical ones. The problem, however, might be more complicated 
because of the possible smaller crystallinity of the spherical particles 
compared to cubes^^ '*^. It is known that surface anisotropy is linked to 
well-defined lattice planes^^, being the spherical entities formed by 
different nanofacets while the cubic particles have fairly flat planes. 

The temperature dependence of zero-field-cooled and field cooled 
(ZFC-FC) magnetization curves establishes further differences between 
the nanomagnetic features of the isolated particles and the collective 
system due to the variations in the magnetostatic interactions^^. For 
example, FC curves (see Fig. S6) reach a plateau in contrast to the 
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Figure 3 | SQUID, (a) Magnetic hysteresis loops recorded at 300 K for 
20 nm and 40 nm square nanoparticles. Inset shows the low field region of 
the hysteresis loop, (b) Magnetic hysteresis loops recorded at 5 K. 

random non-interacting particle system, which points out to a strong 
interaction between particles. The details and the mechanisms behind 
the power loss of particles in such a collective intricate behavior are not 
sufficiently understood yet^^'^^. Indeed we find significant changes in 
the hyperthermia response depending on size, the shape of nanopar- 
ticles and concentration, which do not follow the linear response the- 
ory^. Figure 4 presents the influence of particle concentration on the 
heat production for 20 nm and 40 nm iron oxide nanocubes in aque- 
ous dimethyl sulfoxide (DMSO) solutions. Obviously, the available AC 
field-intensities probe only losses from minor hysteresis loops for the 
40 nm sample (see inset in Fig. 3a), thus explaining the smaller SAR 
values compared to the 20 nm case. We observe also that SAR 
decreases with increasing concentration in both cases and that, as 
expected, it increases in a non-linear way with the maximum applied 
field for all the studied samples. Next, the heating efficiency of 20 nm 
iron oxide nanocubes was compared to spherical nanoparticles, con- 
sidering that both systems are of the same size and crystal structure 
differing only in their shape. Figure 5a shows that the SAR measured 
under the same experimental conditions is about 20% superior for the 
nanocubes, despite similar room temperature soft magnetic features 
are also observed for the spherical nanoparticles (see Supplementary 
Fig. S5 and S6). 

To provide insights into the experimental parameters that dictate 
the pronounced dependence of SAR on concentration and particle 
shape, we performed Monte Carlo simulations (see Fig. 5b). Our 
analysis is based on the macrospin approximation that considers 
each particle as a classical Heisenberg spin with an effective aniso- 
tropy, which can differ from the bulk first-order anisotropy cons- 
tant since it includes surface anisotropy contributions^°'^°'^\ Some 



I 
< 



500 



300 



100 



500 



300 



100 



(a) 



0.5 mg/mL 
1.0 mg/mL 
2.0 mg/mL 





15 20 25 30 

Magnetic field (mT) 

Figure 4 | Hyperthermia. Specific Absorption Rate (SAR) for (a) 40 nm 
and (b) 20 nm iron oxide nanocubes extrapolated from experimental 
thermal response curves at different maximum applied magnetic fields 
(Figure S7 within the Supplementary Information). Data are expressed as 
the mean of 3 measurements ± the standard error of the mean. Inset shows 
suspension stability after measurement. 

representative simulated hysteresis loops are shown also in the 
Supplementary Information (see Fig. S9). Calculations endorse the 
statements of the previous paragraphs: a widely negative infiuence of 
dipole-dipole interactions on the heating power for an increasing 
particle concentration, and the greater hysteresis area for the nano- 
cubes. In view of these observations, it is unambiguous that cubic 
particles exhibit higher anisotropy energy values than that for 
spheres. Although pioneering reports have revealed the imperative 
contribution of surface anisotropy to the magnetic properties of fine 
cubic particles, its role is still an issue of controversy. In order to gain 
some insight on the peculiarities of the experimental magnetic beha- 
vior presented in the previous section, we have conducted a MC 
simulation study of individual spherical and cubic nanoparticles in 
which the magnetic ions forming the particle are considered at the 
atomistic level of detail. We observe that the area of the hysteresis 
loop of the cubic particle is bigger than that of the spherical one (see 
Fig. 6), implying also a higher SAR as observed experimentally. 
Notice that the difference in areas stems from qualitative loop shape 
differences around the coercive and closure field points that can be 
traced back to the different reversal processes of the surface spins (see 
Supplementary Fig. Sll and S12). 

Discussion 

A nice survey of literature, and comparison of experimental and 
calculated SAR values for several nanoparticles types as a function 
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Figure 5 | Comparison of Experimental and Computational Results. 

SAR values for two nanop article solutions of similar concentration 
(0.5 mg/mL) and size volume but different shape indicating enhancement 
of SAR values for the 20 nm square nanoparticles. (a) Experimental results, 
(b) MC simulations for the macrospin model with dipolar interactions at 
300K. 

of the magnetic field frequency (range of 300 kHz- 1.1 MHz and 
with an amplitude up to 31 mT) can be seen in Ref. 10 and 17. 
While commonly available iron oxide ferrofluids show SAR about 
100 W/g^°, in a few special cases experimental values in excess of 
500 W/g^^, and up to about 1 kW/g were found for magneto- 
somes^^'^^. Attention should be paid to the fact that the values are 
usually reported in Watts per iron gram, accordingly a factor about 
1.5 will scale up our results for magnetite, thus concluding that our 
values are among the highest recorded for iron oxide nanocrystals. 
Naively, magnetic losses scale approximately with MsHc oc K given 
that He ~ 0.5Ha within the Stoner-Wohlfarth model, being Ha = 
2 K/Ms the characteristic anisotropy field of the particles (see 
Supplementary Information). Thus, the increased anisotropy in the 
case of cubic particles explains the higher SAR compared to the 
spherical ones. 

The results reported in the present study also reveal that, in 
addition to the surface magnetic anisotropy energy, particle concen- 
tration plays a crucial role in tailoring the heating efficiency inde- 
pendently of the particle geometry, as demonstrated in Fig. 4. 
Reasonably, an increase in concentration corresponds to a decrease 
in the mean interparticle separation and gives rise to a notable 
increase in the dipolar interparticle interactions. The role that dipolar 
interactions might have in SAR is not completely understood at 
present^^ and recent experimental studies have reported either an 
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Figure 6 | Magnetization simulation. Hysteresis loops for a spherical (red 
circles, diameter 20 nm) and a cubic particle (blue squares, side 20 nm) 
obtained from MC simulations of an atomistic spin model of maghemite at 
low temperature. In both, uniaxial anisotropy at the core and Neel surface 
anisotropy have been considered. Snapshots show the spin configurations 
in the remanent state. Spins have been colored according to their 
projection into the magnetic field direction (z axis) from red ( + 1) to blue 
(-1). 

increase or decrease of SAR with interactions^^'^^'^^"^^. Overall, results 
suggest a widely negative influence of dipole-dipole interactions on 
the heating power of nanoparticles, still a clear correlation between 
the magnetic anisotropy of the particles and their hysteresis loop is 
observed and may lead to an increase in magnetic hyperthermia 
efficiency^^. 

Additionally, the different geometrical arrangement of nanopar- 
ticles in suspensions may play some role in explaining the increase of 
SAR for nanocubes compared to spherical particles. For instance, our 
observation of nanocubes chain formation even in the absence of a 
magnetic field (Fig. 2) suggests the existence of strongly anisotropic 
dipolar forces mediating nanoparticle attachments^. We are not aim- 
ing to relate chain formation to concentration increase, but we argue 
that, since nanocubes are more prone to chaining, they would display 
higher SAR values compared to the spherical particles. In this regard, 
a rigorous computation of the magnetostatic energy of chains 
of identical, uniformly magnetized particles of arbitrary shape 
can be found in Phatak et aV^-, it is shown that the face-to-face 
configuration in nanocubes allows for significantly more favorable 
chain formation than the case for spherical entities. Obviously, a 
transient long nanoparticle chain could also form by the application 
of an alternating magnetic field^^, still the self-assembly of colloidal 
crystals into ordered superstructures depends critically on the shape 
(and size) of the nanocrystal building blocks^^. Note that, interest- 
ingly, straight chains made of cubes will have an axial magnetization 
state as the lowest in energy, while it is unlikely that a chain of 20 nm 
iron- oxide nanospheres will be observed experimentally^^. It is 
worth mentioning the hydrodynamic diameters measured by dyna- 
mic light scattering (see supplemental Fig. S3). Results suggest that 
20 nm cuboids self- aggregate in aqueous solution, and the average 
number of particles per chain is roughly 10, confirming the forma- 
tion of agglomerates in these samples while the nanospheres stay 
dispersed. 

Thus, we speculated that the capacity of cuboids to self-assemble 
spontaneously can be used to tailor the heating capabilities. To fur- 
ther cross check this view, we have computed the evolution of the 
hysteresis as a function of the number of particles N within a chain. 
We will focus now on the influence of interparticle dipole interac- 
tions in a chain arrangement as displayed in Fig. 7. Our simulations 
predict that the area of hysteresis loop increases (and therefore of the 
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Figure 7 | Chain-like assemblies. Computed hysteresis loops for arrays of 
nanoparticles of different length (values of N indicated), and the limit case 
of a single particle, where Ha = 2 K/Ms is the anisotropy field of the 
particles. The temperature was introduced in terms of the anisotropy 
energy barrier of the particles, t = kBT/2KV = 0.001. The inset shows the 
magnetic response on increasing the temperature for the N = 1 and 
10 cases, illustrating the higher thermal stability of the chains. The 
schematic picture shows particles within a chain possessing easy axes 
contained within an angle of n/4 (cubes are only for illustration purposes, 
since spherical nanoparticles could also form chains). 

SAR) with the length of the chain. Further, the thermal stability 
gained by creating arrays is an advantage when exploiting hysteresis 
losses (inset Fig. 7). These observations indicate a promising way to 
increase the hyperthermia performance by assembling cubic particle 
in elongated chains. This finding is in remarkably good agreement 
with the results observed for magnetosomes by Alphandery et al.^^ 
Disorientation of the assembly would lead to a considerable decrease 
in the hysteresis loop area and to very different heating properties. 
Example of the later is the decrease on the heating efficiency of 
separate magnetosomes compared with those of magnetosomes 
arranged in chains^^'^^. Independently of the formation of chains 
and its positive effects on SAR, a concentration increase may lead 
to the coalescence of chains leading to cluster formation and a 
decrease in SAR. However, we stress that even in this case, the higher 
SAR of cubes can be accounted for their increased surface anisotropy 
with respect to spheres. 

To conclude, our data and analysis indicate that ferrimagnetic 
nanocubes with an edge length about 20 nm exhibit superior mag- 
netic heating efficiency compared to spherical particles of similar 
sizes. The oriented attachment of magnetic nanoparticles biomi- 
micking magnetostatic bacteria, and the beneficial role of surface 
anisotropy, are recognized as important mechanisms for the 
development of magnetic hyperthermia for cancer treatment. We 
foresee such quantification of nanoparticles interaction and under- 
standings of the magnetization reversal are also important for the 
design of magnetic nanomaterials for other biomedical applications. 

Methods 

Synthesis. Iron oxide nanocubes were prepared by heating a solution of iron(III) 
acetylacetonate (Fe(acac)3), decanoic acid and dibenzylether. This method takes the 
advantage of the discernible separation of nucleation and growth stages caused by the 
intermediate formation of iron(III) decanoate complex as discussed in detail in a 
previous publication*^^. Briefly, size can be tuned over a wide range ( 1 5 nm to 1 80 nm) 
by choosing the suitable heating rate. Namely, for the preparation of 20 nm 
nanocubes (edge dimension), 0.353 g (1 mmol) of Fe(acac)3 was mixed with 0.688 g 
(4 mmol) of decanoic acid in 25 mL of dibenzyl ether. After a short vacuum step at 
60°C (30 minutes), the solution temperature was first raised up to 200°C with a 
constant rate of 2.6°C/min, and kept at this temperature for 2 h under an argon flow 



and vigorous stirring. In a second step, the solution was heated to reflux temperature 
with a heating rate of 6°C/min. After 1 h the solution was cooled down and acetone 
was added. Nanoparticles were then collected by centrifugation at 8000 rpm and 
redispersed in chloroform. This procedure was repeated at least two times in order to 
get rid of the excess of surfactant. Nanocubes of 40 nm were synthesized by 
decreasing the heating rate, during the second step, down to 1.7°C/min. In the 
synthetic procedure for the production of 20 nm (diameter) spherical particles, 

3 mmol Fe(CO)5 were added at 100°C in 10 mL dioctylether in the presence of 

12 mmol oleic acid. The mixture was left to reflux for 3 h at 290°C and then cooled to 
room temperature. Ethanol was added to yield a black precipitate, which was then 
separated by centrifugation. The supernatant was discarded and the particles were 
redispersed in hexane. Given that the resulting magnetic nanoparticles were 
hydrophobic, the powders were further dissolved in a mixture of dimethyl sulfoxide 
(DMSO) and water (1:1). DMSO is a naturally derived, inexpensive, non-toxic 
solvent and pharmaceutical agent that has been demonstrated to be a well-tolerated 
excitatory modulator in the management of cancer pain^^. Furthermore, DMSO has 
chemical properties which facilitate its absorption into and distribution throughout 
biological systems by all routes of administration, thus biocompatibility in future 
applications is guaranteed. 

Characterization. Iron oxide particles were characterized by transmission electron 
microscopy (TEM) using a JEOL JEM-2100 (200 keV) for high resolution (HR)TEM, 
and by a field-emission gun scanning transmission electron microscope (FEG)TEM 
FEI Tecnai F20-G2 (200 kV) equipped with a high-angle annular dark-field detector 
STEM-HAADF for three-dimensional (3D) electron tomography. The samples were 
prepared by dropping on solution of nanoparticles onto a carbon coated copper grid. 
The 3D structure of nanoparticles, as well as that of the assembled chain, is 
reconstructed from a tilt series (range of —64° < a < 64° at an increment of 

4 degrees) of 2D projections, using the simultaneous iterative reconstruction 
technique (SIRT). In order to determine the apparent hydrodynamic diameter, 
dynamic light scattering (DLS) analysis was carried using Malvern Instruments 
Hydro 2000MU accessory. Thermogravimetric analysis (TGA) was carried out with a 
TGA-SDTA 851e/SF/l 100 Mettier Toledo device up to 700°C, by heating the sample 
under an argon/nitrogen flow at a heating rate of 10°C/min. The residual weight 
accounts for the mass of iron-oxide nanoparticles in the ferrofluid. Quasi-static 
magnetic characterization was carried out in a superconducting quantum 
interference device (SQUID) Quantum Design MPMS XL-7T magnetometer. M(H) 
hysteresis loops were measured at different temperature (5K and 300K) by applying 
field up to 5 Tesla. Magnetization zero field cooled (ZFC) measurements were 
performed upon warming with an applied magnetic field of 5 mT after cooling the 
samples in zero applied field. The field cooled (EC) curves were obtained by 
measuring at stepwise-decreasing temperatures in the same small applied field. AC 
magnetic hyperthermia experiments were performed using a 23 mm diameter three- 
turn induction coil powered by a 4.5 kW AC field generator. While frequency was 
kept constant at 765 kHz, the amplitude of the applied magnetic field was tuned from 
1 5 up to 30 mT. The amplitude of the alternating magnetic field was estimated using a 
pick-up coil connected to an oscilloscope. Temperature was monitored by using an 
OpSens PicoM GaAs-based fiber optic probe immersed in a test tube containing 2 mL 
of solution. Three different iron-oxide concentrations (0.5, 1.0 and 2.0 mg/ml) were 
exposed to the alternating magnetic field for 900 seconds. Specific absorption rate 
(SAR) values were estimated by subtracting the solvent background signal and the 
heat losses to the environment. Further details on the hyperthermia capabilities of 
these particles under safe clinical conditions (Hjnax/< 5 X 10**A/ms), and 
intracellular uptake trials, can be found in a recent publication^ \ 

Computational details. The simulations were performed both at the atomic level and 
under the so-called macrospin approximation, in order to investigate the single- 
particle properties in relation to shape and anisotropy, and the role of magnetic 
dipolar interactions, respectively. First, we used the Monte Carlo (MC) method with 
the standard Metropolis algorithm. The physical model employed for our numerical 
simulations considers a perfectly monodisperse system of single-domain magnetic 
particles with an effective uniaxial magnetic anisotropy. Though, it should be pointed 
out that the SAR may be also influenced by the particle size distribution (as recently 
considered by Hergt et aV"^ and Carrey et al. ^°), which, however, is beyond the scope of 
the present paper. The spatial distribution of the particles resembles a frozen 
ferrofluid without aggregations, the positions of the particles are kept fixed and the 
easy axes are chosen randomly. The energy model is the same as in Ref. 48 so that the 
energy of each particle in this ideal scenario is regarded to have three main sources: 
anisotropy, Zeeman and dipolar interaction. In our simulations, we reproduce M(H) 
curves at different sample concentrations and for different values of maximum 
applied field Hmax- To fit the simulations with the experimental procedure, we applied 
magnetic fields of the same amplitude as were applied experimentally. The model 
shows that SAR arrives at the summit for large field amplitudes, for values about the 
characteristic anisotropy field of the particles Ha = 2 K/Ms. Next, we have 
considered an atomistic model with Heisenberg classical spins placed at the magnetic 
ion lattice sites of the real maghemite structure (see Supplementary Information for 
details about the model and simulations techniques), and we have computed the low 
temperature hysteresis loops of two individual particles with spherical and cubic 
shapes having the same size as the ones experimentally studied. Finally, in order to 
investigate the degree of magnetostatic coupling between nanoparticles, we studied a 
series of chains with N particles. The anisotropy axis of the particles was arbitrary 
taken within the chain axis, equally distributed apart from the median position in 
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order to account for deviations in the positions of the crystals from the chain axis, thus 
resembling the results obtained by electron holography on magnetosomes*^^. 
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